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ABSTRACT 
We    derive    the    probability   distribution    P{Xn  =  x}   for   a   random    walk    satisfying 
X0=0;   P{Xn  +  l  =  (2X„+1)  modp}  =  P{Xn  +  1  =  (2XB-1)  mod  p}  =  y  for  n  >  0 
and  p  odd  integer. 

Chung  et.  al.  [1]  present  a  scheme  for  combining,  or  shuffling,  random  number  gen- 
erators of  the  form 

ATn  +  1  ■  aX„  +  b(mod  p),     p  odd;      Xq  —  0. 

Their  scheme  consists  of  choosing  b  from  another  truly  random  generator,  in  particular 
they  deal  with  the  following  three  cases: 

I:        a  =  1,      b  =  0,  1,  -1       each       w.p.  — 

II:       a  =  2,      b  =  0,  1,  -1      each      w.p.  — 

III:       a  =  2,      b  =  1,  -1  each      w.p.  — 

In  each  case,  they  furnish  the  number  of  steps  which  suffice  to  guarantee  randomness, 
after  obtaining  the  variation  distance  (see  (1.6))  as  a  measure  of  closeness  to  uniformity. 

The  purpose  of  this  note  is  to  give  a  simple  derivation  of  the  distribution 

P{Xn  =  x}     x  =  0,...,p-\ 
for  case  III,  from  which  the  proof  of  Theorem  5  in  [1]  is  easily  obtained. 

Theorem: 

Let  X0,X1PC2,  •••  be  a  sequence  of  random  variables  satisfying 

(1.1) 

P  {Xn  +  1  =  2X„  +  \  modp}  =  P{Xn  +  1  =  2X„-\  modp}  =  y 

for  nSO  and  p  odd  integer. 
Let 

2"  =  ap  +  p  +  1      0<p  <  p-\ 


Then  for  x  =  0,...,  p  —  1 


P{Xn  =  x}  = 


a  +  1 


2" 
a 
2" 


ifx=(p-2;)modp       ;  =  0 p 


Proof: 

For  n^l  define 


(1.2) 


X„  =  2""1  o-j  +  2""2ct2  +  ...+  cr„    morfp 
where  a,-,    j=1,...,/i  are  i.i.d  and 

P{a.=  +i}  =  P{a.=  _1}=  1  . 
Then,  the  sequence  (1.2)  satisfies  (1.1)  and 


(1.3) 


^    "         =    U  E  e    p  \=    n 


i  u^r2  .    i  -ik~^2 


f  sin2"  2tt£/p 


n  cos  it 2' 

/=o  p 


2"sin  lirk/p 


p       +  4. 


//    it    *    0 


if  k    =    0 


) 


But 


1/>-l      -Or x       [    i* X.\ 


(1.4) 


p 


p-    -A 


1+2  * 
Jt=i 


ik^-X. 
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2Trk                         2Ttk 
For  2"  =  ap  +  P+1     0<(3<p-l      sin  2"    =  sin((5  +  l) hence,  using  (1.3) 


and  (1.4)  we  obtain 


P{Xn  =  x}=± 


,      ..2v    sin(B+l) 

1    p  —  1      —  uc X  v  "  ' 


2tt£  1 


1  +  TTS   e 

£    k=\ 


sin 


2ir* 


2"p 


2"  -  (p  +  1)  +   X   B        P    


2-nk 


2np 


2-nik  .„     _  ,      . 
P-l     — —  (P-2j-*) 

ap  +  2    2  e  p 

j=0    *  =  0 


or 


(1.5) 


P{Xn  =  x}  =    —  jap    +    J   P8*,P~27  -Omodp  \ 

2  P    I  j=0  J 


a+1 


2" 


if  x=($-2j)  mod  p        ;  =  0 0 


—  otherwise. 

2" 


for  jr  =  0,...,/J  — 1. 

The  variation  distance  is  defined  in  [1]  as 


Q.E.D. 


(1.6) 


Z    x  =  0 


P{Xw=x}- 


o.          1     _.   a  + 1         .  ,.  <•  /i  c\  ■ 

since  —  <  —  <  an  immediate  consequence  of  (1.5)  is 

2"        P  2" 
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z    x=0 


which  is  Theorem  5  in  [1]. 


P{Xn=x}-j 


=    (P+D[P-(P  +  1)] 
2"p 


1.  Acknowledgement 

I  would  like  to  thank  Malvin  Kalos  for  introducing  me  to  the  problem  that  motivated 
this  work. 

2.  References 

[1]     Chung,  F.R.K.,  P.  Diaconis  and  R.L.  Graham.  "Random  Walks  Arising  in  Random 
Number  Generation",  Annals  of  Probability. 


Ultracomputer  Note  113 


Page  4 


NYU   UCN-113  C.l 

Percus,  Ora  E 
A  note  on  random  number 
generator  of  Chung  et  al. 


NYU   UCN-113  c.l 

Percus,  Ora  E 
A  note  on  random  number 
generator  of  Chung  et  al. 


This  hook  may  be  kept         jyjy     q    j     jQQ-j 

FOURTEEN    DAYS 

A  6ne  will  be  charged  for  each  day  the  book  is  kept  overtime. 

GAYLORD   142 

PRINTED   IN   u    S    A 

